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Knowledge of taxis (directed swimming) in the Archaea is currently expanding through identification of 
novel receptors, effectors, and proteins involved in signal transduction to the flagellar motor. Although the 
ability for biological cells to sense and swim toward hydrogen gas has been hypothesized for many years, this 
capacity has yet to be observed and demonstrated. Here we show that the average swimming velocity 
increases in the direction of a source of hydrogen gas for the methanogen, Methanococcus maripaludis using 
a capillary assay with anoxic gas-phase control and time-lapse microscopy. The results indicate that a 
methanogen couples motility to hydrogen concentration sensing and is the first direct observation of 
hydrogeno taxis in any domain of life. Hydrogenotaxis represents a strategy that would impart a competitive 
advantage to motile microorganisms that compete for hydrogen gas and would impact the C, S and 
N cycles. 



Hydrogen gas (H2) is a crucial substrate for methanogens as well as a common source of energy for other 
organisms in both anaerobic and aerobic environments, including acetogens, sulfate- and sulfur- reducers, 
and hydrogen-oxidizers^"^. Biological methane (CH4) production from H2 and carbon dioxide (CO2) 
contributes to greenhouse gas emissions and is possibly one of the oldest microbial metabolisms^'^. Under- 
standing the ecological strategies of methanogens is not only important for our knowledge of early earth processes 
and present-day anaerobic environments, but also for determining potential roles in human health conditions 
{e.g., colon cancer and periodontal disease), where positive correlations have been made with incidence of disease 
and occurrence of methanogens^'^. Methanococcus maripaludis is an anaerobic archaeum that can use H2 or 
formate as electron donor to reduce CO 2 to CH4 and is considered a model mesophilic methanogen. Recently, the 
swimming behavior of M. maripaludis was described^, but chemotactic responses have not been shown. 
Chemo taxis has been demonstrated for Archaea, including methanogens but taxis to hydrogen has not been 
shown for any domain of life. The chemotaxis signal transduction system in Archaea is similar to the well-studied 
system in Bacteria; however, the flagellar switch is different and none of the archaeal flagellar proteins have 
homologs to bacterial flagellar proteins 

Chemotaxis has been the subject of many mathematical models and the majority have concentrated on 
reproducing the population -level observation of migrating bands of high cell concentration in swarm plates 
and capillary experiments^^. Pioneering work in modeling chemotaxis behavior by Keller and Segel in 1971 has 
been the basis of the most common mathematical models^^. In one dimension, with x being the spatial variable, 
the Keller-Segel model can be described as a flux, J, such that 



J = 



db ^ ds 



dx 



dx 



(1) 



where |i is the cell diffusion coefficient that takes random, non- directed, movement of cells into account, b is the 
microbial population density, s is the attractant concentration and x(s) is the non-constant chemotactic coef- 
ficient. The population flux, J, can be differentiated to yield the more common form 
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and the average cell swimming velocity, v, is calculated by dividing the flux by the population density, or v = J/b. 
Lapidus and Schiller^^ proposed a form of x(s) such that 
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where % is the constant chemotactic coefficient and is the receptor- 
Ugand binding dissociation constant. The Lapidus- Schiller x(s) term, 
and variations thereof, have been used widely to describe chemotaxis 
in bacteria^^. Most work on archaeal chemotaxis has been performed 
with Halobacterium salinarum; however, attempts to mathemat- 
ically describe population flux in archaea have not been published. 
In addition, % and kj have not been determined for any archaea^^'^°"^^. 

The goal of the present study was to subject M. maripaludis cells to 
a H2 concentration gradient and compare swimming behavior to 
model predictions. Attractant (H2) transport was modeled by 
Fickian diffusion and consumption by the population was modeled 
in the Michaelis-Menten form such that 

dt~ dx^ (km+s) ^ ^ 

where D is the diffusion coefficient for H2, r^ax is the maximum 
consumption rate, and k^ is the half- saturation constant. 

A modified capillary assay was used in which cells were loaded into 
a gas-tight anaerobic capillary under an anoxic atmosphere, and a 
valve allowed controlled exposure to a H2 source (see Supplementary 
Fig. SI and S2 online). In previous experiments, capillaries that con- 
tained a dissolved chemoattractant were immerged into a cell sus- 
pension and cells entered the capillary in the presence of a 
chemoattractant over a known incubation time. The population 
increase in the capillary relative to a control was then quantified 
through cell enumeration^^. For this study, a novel method was 
developed to directly observe microscopic swimming behavior of 
anaerobic cells inside a capillary during exposure to a gas. M. mar- 
ipaludis cells were tracked for direction and velocity changes upon 
exposure to H2 or an Argon (Ar) control to quantify the average 
population-wide response. Cell movement was measured using 
time-lapse confocal laser scanning microscopy (CLSM) in the center 
of the capillary (0.5 cm from the gas phase; as in Supplementary Fig. 
SI online). A reaction -diffusion model predicted that within 10-15 
minutes H2 concentrations would reach the threshold of 2.5-23 jiM, 
at which hydrogenotrophic methanogens have been shown to use H2 
in pure culture studies^^"^^ (Figure 1). 

Results 

When M. maripaludis cells were exposed to H2, the average swim- 
ming velocity exhibited significant bias toward H2 within ten min- 
utes (Figure 2A). Biased random walk was not observed when cells 
were exposed to an Ar gradient (Figure 2B), nor was swimming 
velocity affected normal to the H2 or Ar gradients (Figures 2C and 
D). A strong chemotactic response was observed when cells were first 
starved of H2 for 4-5 hours, while cultures that had not been starved 
did not show an increase in biased swimming (see Supplementary 
Fig. S3 online). Similar results were observed in Rhodobacter sphaer- 
oides and Sinorhizobium meliloti when well-fed cells exhibited 
weaker chemotactic responses than starved cells toward attractants 
such as organic acids and sugars. It was suggested that this was 
because chemosensory pathways in these organisms were dependent 
on metabolic state and/or transport^^, as could be the case for H2 
sensing in M. maripaludis. 

The starvation period (4-5 h) that was required to induce a tract- 
able response to H2 did not cause loss of motility, as observed during 
starvation of H. salinarum^^. The average swimming speed of starved 
M. maripaludis cells was 2.1 ± 0.05 jim s"^ before exposure to H2 
(Figure 3A) and increased to 3.1 ± 0.02 jim s"^ after exposure; equal 
to non- starved cells. Maximum swimming speed averages (calcu- 
lated by averaging the single maximum speed from each time point) 
were highly variable between time points, and the average maximum 
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Figure 1 | The predicted hydrogen concentration over time at the 
observation point 0.5 cm from the gas phase over the course of the 
experiment shown with (A) linear axes and (B) log-linear axes. The 2.5- 
23 |aM threshold, at which hydrogenotrophic methanogens have been 
shown to use H2 in pure culture studies^^"^^, is reached at approximately 10 
minutes. 

was higher for non-starved cells (9.3 jim s"^ versus 8.5 jim s"^ for 
starved cells + H2) and lowest for starved cells before H2 exposure 
(7.0 |im s~^) (Figure 3B). The highest maximum observed speed was 
91 jim s~^ for H2 starved cells after exposure to H2 (Figure 3B inner 
boxes), approximately twice the previously observed maximum 
speed of 45 jim s"^ for M. maripaludis^. 

A Keller-Segel model was applied with boundary conditions spe- 
cific to our study in an effort to quantify the observed chemotactic 
response to H2 in context with other known chemotactic responses. 
Currently, K-S model parameter values do not exist for any organism 
in Archaea or for any other gas besides O2; therefore, a broad range of 
kd values were used in the model. Three unknowns, namely ligand- 
receptor dissociation constant (k^), chemotactic coefficient (x), and 
random cell diffusion coefficient (|i), were independently fitted by 
varying one unknown across the range of published literature values, 
while keeping the other two variables constant at the average lite- 
rature values (see Supplementary Table SI online). The best fit 
for average swimming velocity was obtained with a k^ value of 
0.70 mM, with 0.30 and 2.30 mM corresponding to the 95% confid- 
ence interval of the data, while % and \i were kept at average published 
values (Figure 4A). This range for k^ is similar to that observed for 
Escherichia coli AW405 to a-methyl aspartate^^. It is; however, quite 
different from values reported for Bacillus subtilis receptor affinity to 
O2 (0.0015 and 0.075 mM for the high and low affinity of the recep- 
tor, respectively), which is the only previously reported k^ for a gas^°. 

The model could also be fitted to the velocity curve by varying % 
and keeping k^ and ji constant at average literature values 
(Figure 4B). The average % of 9.3 X 10"^ cm^ s~^ used to fit the 
experimental results is higher than the literature range 7.20 X 10"^ 
to 1.24 X 10"^ cm^ s"^ (see Supplementary Table SI online). The 
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Figure 2 | Average cell velocity (black) with 95% confidence intervals (gray) before and after opening valve to gas phase for (A) H2 and (B) Ar control. 

Positive y-axis values indicate movement toward the gas phase, negative y-axis values indicate movement away from the gas phase. Average cell velocity 
(C) and (D) normal to concentration gradient (n = 5 for H2 and n = 3 for Ar). 



experimental results were best fitted to the model that assumed the 
average published value for and then by varying (Figure 4A). 
The model could not be fitted to the experimental data by varying ji 
according to published values (Figure 4C, note the difference in y- 
axis scale); however, changing ji allowed for a change in shape of the 
velocity curve. 

The response time predicted by the model does not fit the rapid 
response observed in the experiment and this is likely due to an 
inaccurate prediction of mass transport inside the capillary (for both 
H2 and cells). Both random cell diffusion (ji) and the diffusion of 
hydrogen would require correction from predicted values if mass 
transport was inaccurately represented. First, considering |i indepen- 
dently, it is likely that response time would be faster for cells with a 
larger ji (Figure 4C). This is based on the proportionality of |i to 
swimming speed and run time, and inverse proportionality to one 



minus the cosine of the turn angle Although not measured 
directly here, M. maripaludis has been shown to have relatively long 
runs and small changes in direction or turn angles^. Turn angles 
between 20 and 45° would result in the largest value of |i for a given 
swimming speed and run time, so it is reasonable to assume |i would 
be higher for this type of swimming. The observed rapid response can 
be explained by an inaccurate prediction of H2 mass flux into the 
hquid domain. If H2 were to reach the observation point faster than 
predicted by the model, then one would predict a proportionally 
faster response. This was an entirely static system on a visually 
observable scale but it is possible that unpredictable micro-scale 
convective forces increased H2 mass transport. The collective, 
directed motion of swimming cells may have induced convective 
flow inside the capillary. Previous work estimates that a force of 
0.5 pN is exerted by a cell swimming at 25 jim s~^ While 




Figure 3 | Swimming speed of M. maripaludis cells when starved for 5 hours before H2 exposure ( — H2 and n = 4,725), after H2 exposure ( + H2 and n = 
20,189), and when not starved (No starve, n = 45,001). Values represent the mean and error bars represent 95% confidence intervals. (A) Average 
swimming speed. The difference between -H2 and both +H2 and no starvation is significant (p < 0.05). (B) Average maximum swimming speeds are 
represented by bars, and absolute maximum values inside each bar represent the highest observed speed over all time points for the condition. Differences 
between all conditions are significant (p < 0.05) and p values were calculated with two-tailed t-tests. 
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Figure 4 | Model fitting results (smooth black lines) overlaid on experimental results with 95% confidence intervals (grey lines) (A) using average 
literature values for % and \i and varying over the range shown or (B) using average literature values for and \i while varying % over the range shown 
(C) or \i varied across the range of literature values. 



chemotaxis -induced convective flow is difficult to demonstrate con- 
clusively, increased mass transport in the capillary could affect mass 
transport predictions. 

A diffusion coefficient correction term, 5, was introduced as a tool 
to investigate the effect of increased mass transport on the predicted 
response time. The correction was applied to the H2 diffusion coef- 
ficient and the random cell diffusion coefficient (|i) simultaneously, 
making the assumption that micro -scale convection was the primary 
source of error. Any convection would affect the random diffusion of 
cells and the diffusion of H2 equally. Increased mass transport in the 
Hquid (simulating convection) caused an earlier response in all con- 
ditions that were investigated (see Supplementary Fig. S4 online). 
The value of 5 is likely to be on the order of 1 X 10"^ to 1 X 
10"^ cm^ s"^ based on curve shape and response time as compared 
to the experimental velocity data. 

Another possible explanation for the more rapid experimental 
response is the presence of more than one type of H2 receptor with 
varying affinities. This would allow M. maripaludis to respond to H2 
across a wider range of concentrations. The predicted kj value 
0.7 mM is high, and the presence of a high affinity receptor in M. 
maripaludis is likely. B. subtilis O2 receptors with two distinct affin- 
ities and binding components have been shown^°. The presented data 



demonstrate a chemotactic response to H2, but more focused physio- 
logical studies and subsequent model refinement are needed to better 
represent this poorly understood phenomenon. 

In our capillary assay, cells were not observed to accumulate in 
bands despite the chemotactic response observed. Typical chemo- 
tactic bands observed in capillary and swarm plate assays result from 
cell accumulation as a net result of a biased random walk^^. In the 
capillary assay used, convection inside the capillary may have pre- 
vented cell accumulation. It is also possible that there are factors 
involved in chemotactic band formation that are unique between 
organisms or factors specific to hydrogenotaxis {e.g., swimming 
mode, adaption response, quorum sensing). 

In our work, banding-like behavior was observed on a larger scale 
when M. maripaludis batch cultures were grown statically with H2. 
Under these conditions a pellicle formed at the gas-liquid interface 
(Figure 5 A, B). This is a similar observation to the one made by 
Beijernck in 1893, where aero tactic cells were observed swimming 
toward the meniscus of a test tube^^. When cultures of M. maripa- 
ludis were grown statically with the soluble electron donor formate, 
cells grew throughout the liquid medium (not in a pellicle) and had 
less cell-associated carbohydrate than H2-grown pellicle cultures 
(see Supplementary Fig. S5 online). These results suggest that 
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Figure 5 | (A) and (B) Field emission scanning electron micrographs of M. maripaludis pellicle with extracellular material and cells (indicated by 
arrows in (A)). 



extra-polymeric substance (EPS) production is important for the 
formation of a pellicle in response to a H2 gradient and may require 
a time-scale longer than used in the capillary experiment. 

Discussion 

The hydrogenotrophic methanogen, M. maripaludis, displays che- 
motactic behavior toward H2 gas (hydrogenotaxis). Changes in pH 
were not observed in medium with or without H2 under the tested 
growth conditions and times, and these results indicate that the 
observed taxis was in response to H2 and not H^. Although the 
commonly used Keller-Segel model was not able to predict the exact 
response with previously published parameter ranges (values), para- 
meter adjustment allowed the model to replicate the trend on a scale 
similar to other organisms and chemoattractants. The ability to move 
toward higher concentrations of H2 could incur an advantage to 
methanogens that are otherwise outcompeted by those that are able 
to use H2 at lower concentrations and/or utilize terminal electron 
acceptors that are more energetically favorable. The demonstrated 
chemotactic response would also allow cells to maintain desirable 
localization with respect to the major energy source as well as allow 
for proximity to H2 -producers in mixed communities. Thus, hydro- 
genotaxis could play a crucial role in the establishment and 
maintenance of microbial interactions at the population- and com- 
munity-level. The observed hydrogenotaxis could represent a wide- 
spread eco-physiological strategy of methanogens and other 
hydrogen-utilizing microbes that are important to processes such 
as bioremediation and overall carbon cycling^'^l CH4 is a potent 
greenhouse gas that has an estimated global warming potential 
(GWP) 25-40% higher than CO2 per molecule'^ The three largest 
contributions of CH4 to atmospheric flux as of 2010 (wetlands, 
ruminant emissions, and rice cultivation) are the net result of the 
activity of anaerobic communities dominated by the exchange of 
p^^35,36 ^Yie best of our knowledge this is the first direct observation 
of hydrogenotaxis in any domain of life. 

Methods 

Culturing conditions. Methanococcus maripaludis strain S2 was grown in Balch 
tubes or serum bottles fitted with black butyl stoppers (Geo -Microbial Technologies 
Inc., Ochelata, OK) and aluminum crimp seals. Methanococcus Culture Medium 
(MCC) was prepared under a stream of anoxic 80% N2, 20% CO2 and contains per 
liter 0.33 g KCl, 2.7 g, MgCl2»6H20, 3.5 g MgS04-7H20, 0.14 g CaCl2-2H20, 0.5 g 
NH4CI, 5 g NaHCOs, 22 g NaCl, 0.14 g K2HPO4, 5 mL FeS04 solution (0.19 g 
FeSO4-7H2O/100 mL of 10 mM HCl), 1 mL trace metal solution (per 100 mL; 2.1 g 
Na3Citrate-2H20, adjust pH to 6.5, then 0.45 g MnS04-H20, 0.1 g CoCl2-6H20, 
0.1 gZnS04-7H20, 0.01 g CuS04-5H20, 0.01 g A1K(S04)2, 0.01 gH3BO4,0.1 g 
Na2Mo04-2H20, 0.025 g NiCl2-6H20, 0.2 g Na2Se03, 0.01 g VCI3, 0.0033 g 
Na2W04*2H20) 10 mL of vitamin solution (per liter; 2 mg biotin, 2 mg folic acid, 
10 mg pyridoxine HCl, 5 mg thiamine HCl, 5 mg riboflavin, 5 mg nicotinic acid, 
5 mg DL-calcium pantothenate, 0.1 mg vitamin B12, 5 mg p-aminobenzoinc acid, 
5 mg lipoic acid), 1 mL of Resazurin solution (1 g/L). This solution is boiled under a 
stream of gas before adding 0.5 g cysteine*H20»HCl, then cooled under gas^^. The 
solution was dispensed anaerobically and autoclaved. After inoculation, the 
headspace was displaced and pressurized to 25 PSI with anoxic 80% H2:20% CO2 
through a sterile filter. Modified MCC medium was used for growth experiments with 



formate, where NaCl was reduced to 10.5 g/L; 200 mM Na-formate and 200 mM 3- 
(N-morpholino) propanesulfonic acid (MOPS) were added, and no H2 was added. 

Cultures for taxis experiments were grown at 30°C (with no shaking) from frozen 
glycerol stocks in 40 mL MCC in 125 mL serum bottles to late exponential phase, 
then transferred to 5 mL MCC (10% inoculum) in 18 X 150 mm Balch tubes. These 
working cultures were grown to stationary phase (around 120 h and average cell 
density 69,420 cells/mL). 1 mL of culture was added to 5 mL fresh MCC in a Balch 
tube with no H2 and either incubated for 4-5 h at 30°C under starvation conditions 
(no electron source), or used immediately for non- starvation conditions. Media salts 
were separated from cell suspension prior to taxis experiments by centrifuging the 
inverted Balch tube (100 X g; 30°C for 10 minutes). 

Electron microscopy. Pellicles were removed from tubes with a plastic inoculating 
loop and placed on a round coverslip freshly treated with poly- L- Lysine (1 |ig/mL) 
and fixed in a solution of 2% paraformaldehyde, 2.5% glutaraldehyde and 0.05 M Na- 
cacodylate overnight at room temperature. Coverslip s were rinsed and stepwise 
dehydrated in ethanol, and then critical point dried on a Samdri-795 (tousimis, 
Rockville, MD). Coverslips were mounted on SEM stubs with double-sided carbon 
tape and colloidal silver, and then sputter coated with Iridium for 35 s at 35 mA. 
Images were collected on a Zeiss Supra55VP FE-SEM. 

Capillary assay. Square glass capillary tubes (1.0 mm) with the ends fitted with 
norprene tubing connected to a polypropylene female luer-lock hose barb adapter 
(Cole-Parmer, Vernon Hills, IL) were partially filled with the cell suspension in an 
anaerobic chamber that contained only N2 and CO2. A 5 mL glass gas-tight luer lock 
syringe (SGE, Inc., Austin, TX) was used to transfer the cells to the capillary, and left 
attached to the capillary with the valve closed. A second 5 mL glass gas-tight syringe 
with 100% H2 was attached to the gas side of the capillary (see Supplementary Fig. SI 
online). The entire syringe/capillary assembly was removed from the anaerobic 
chamber and placed in a petri dish water bath on the microscope stage and firmly 
secured with poster putty and tape (see Supplementary Fig. S3 online). 

Microscopic observation of swimming behavior and image analysis. A Leica TCS 
SP5 II upright confocal microscope was enclosed in an incubation chamber and was 
equilibrated to 30° C. High-resolution time lapse images were collected every 0.753 s 
at the center of the capillary, 0.5 cm from the cell suspension/gas phase interface. A 25 
X water-dipping objective was used and a 3 X optical zoom was applied resulting in a 
final field of view of 206.9 X 206.9 [im^ and a pixel size of 0.20 |im. Images were 
initially acquired for 10 minutes with N2/CO2 in the gas phase, then the valve was 
opened to the H2 syringe and images were captured for 40 minutes. The control 
experiments were identical to the above except that 100% Ar was used instead of H2, 
and images were acquired for a minimum of 8 minutes before the valve was opened 
and 38 minutes after. 

Images were manually thresholded and binarized using MetaMorph v. 7.6 
(Molecular Devices, Sunnyvale, CA). Binary images were analyzed using Imaris v. 
7.5.2 (Bitplane, Inc., South Windsor, CT) with a particle-tracking module (Imaris 
Track). 1 s, 5 s, and 10 s filter durations were tested where a given track was only 
analyzed if it was as long or longer than the specified duration. The 5 s filter was used 
for the described analyses, and differences were not observed in overall trends 
between track lengths. 

Chemotaxis model. A one dimensional finite element model was constructed using 
Comsol Multiphysics Version 4.3 a that solves Equations 2 and 4 simultaneously in 
the liquid domain. Diffusion of hydrogen through the gas domains was modeled by 
Maxwell- Stefan equations. All model parameters were corrected for temperature 
(30°C) and salinity (2.65% m/v) of the medium, where possible, and hydrogen 
consumption rate parameters were estimated from literature values. Supplementary 
Table SI online shows all constants used in the model. 

The geometry of the model consisted of one liquid domain and two gas domains 
separated by a valve that opens at t = 0 to start the diffusion of hydrogen into the 
system. The short segment between the valve and the far right boundary is the length 
of the connection between the valve and the main volume of the gas-tight syringe. The 
geometry for each experimental replicate was slightly different so average lengths 
were used. The diffusion of hydrogen through the gas domain of the capillary was 
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expected to be much faster than diffusion through the liquid domain so the precise 
length of the gas domain was not expected to affect the analysis. The diffusion- 
limiting segment between the liquid-gas interface and the point of observation 0.5 cm 
into the liquid remained constant through all replicates and is therefore not an 
average value. Initial conditions are zero hydrogen inside the capillary, 100% 
hydrogen behind the valve in the second gas domain and a constant population 
density in the liquid domain. Supplementary Figure SI online shows the geometry, 
initial conditions and boundary conditions used in the model. 

The model geometry was meshed at a maximum element size of 0.1 mm and time 
step discretization was done with a backward differentiation formula (BDF) method. 
The time steps taken by the solver were allowed to be free with larger steps being taken 
at later time points where gradients are shallower. A relative solution tolerance of 
10"^" was imposed and PARDISO^** was the general direct solver as implemented in 
Comsol. 

Most parameters in the model are well known or able to be calculated from the 
literature with the exception of the chemotaxis constants x> 1^ and k^. A literature 
search was performed for observed chemotaxis values in any organism to establish a 
range of reasonable values. Twelve applicable values were found for twenty-one for 
|i and seven for k^. Maximum, average and minimum values for each are found in 
Table SI. The model was fitted to the average swimming velocity data by indepen- 
dently varying one chemotaxis parameter by trial and error while keeping the other 
two constant at the average literature value found in Supplementary Table SI online. 

Every effort was made to model the diffusion of hydrogen into the liquid domain 
accurately; however, the likelihood of unpredictable factors such as micro- scale 
mixing due to convection still exists. To provide the flexibility to correct for enhanced 
mass transport, a correction term was applied such that 

D'Ha,=DHa,+S (5) 

Where 6 represents the mass transport enhancement beyond what is predicted from 
Fickian diffusion alone and D'naq is the corrected diffusion coefficient used in the 
model for this analysis. Similarly, 6 was also applied to the random cell diffusion 
coefficient such that 

pi' = fi-\-S (6) 

because any correction applied to Dnaq would need to be applied to [l on grounds that 
convection would affect the movement of cells the same as the pass transport of H2. 

Carbohydrate and protein measurements. Protein concentrations were determined 
with the Lowry assay using bovine serum albumin as the standard^''. Hexose sugars 
were measured by the colorimetric cysteine- sulfuric acid method with glucose as the 
standard. Pentose sugars were measured with a colorimetric orcinol-FeCla assay with 
xylose as the standard. A colorimetric carbazole assay was used to measure uronic 
acid concentration with D-galacturonic acid as the standard"^". 
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